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Abstract 

We perform a momentum-space calculation of the reaction pp ppir near thresh- 
old. Direct production and pion rescattering are taken into account. For the latter 
a T-matrix obtained from a microscopic model of the ttN interaction is employed. 
The effects of some approximations commonly used in studies of pp — ► ppir are dis- 
cussed. We find that those approximations overestimate the cross section and also 
modify its energy dependence. We obtain a significant contribution from rescatter- 
ing, but not nearly enough to describe the data. 



1 Introduction 

The high precision data of the reaction pp — > ppir close to pion produc- 
tion threshold from the Indiana University Cyclotron Facility (IUCF) [1,2] 
continue to attract the attention of theorists. Early theoretical studies [3,4] 
followed closely the conventional approach to s-wave meson production which 
goes back to the work of Koltun and Reitan from 1966 [5]. In this model 
two production mechanisms are considered: (i) direct production of a pion, 
depicted in Fig. la, (ii) production of a pion which first scatters off the other 
nucleon before emission, as shown in Fig. lb. It became clear that such a 
model grossly underestimates the empirical cross section [3,4]. On the other 
hand, the predicted energy dependence, which in this model is determined 
essentially by phase-space factors and the pp final-state interaction [1], was 
found to be in agreement with the IUCF data. 

A first quantitative description of the pp — > ppir cross section was achieved 
by Lee and Riska [6]. These authors considered in addition effects from meson- 
exchange currents due to the exchange of heavy (i. e. short ranged) mesons, 
as shown in Fig. lc. The resulting contributions (in particular that of the a 
meson) enhance the cross section near threshold by a factor of 3-5 [6,7] and 
thus eliminate most of the underprediction found in earlier investigations. As 
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a consequence it was argued in ref. [6] that the reaction pp — > ppir could 
provide direct information on the short range part of the NN force. 

However, in a very recent investigation by Hernandez and Oset [8] an al- 
ternative explanation for the missing strength in the 7r° production close to 
threshold is presented. Following the idea of Hachenberg and Pirner [9] the off- 
shell properties of the ttN amplitude are taken into account explicitly when 
evaluating the rescattering diagram. Since the off-shell ttN isoscalar ampli- 
tude is, in general, much larger than its on-shell value, which was used in 
former computations, the contribution from rescattering is now considerably 
enhanced and the resulting pp — > ppir cross section can even be brought 
into agreement with the experimental data. Thus their result suggests that 
the reaction pp — > ppir might be a tool to pin down the off-shell behavior of 
the ttN interaction. 

In the present paper we re-examine critically some of the published results. In 
these publications several approximations to the pion production vertex and 
the kinematics have been tacitly taken over from the very first investigation 
by Koltun and Reitan [5] and we want to shed some light on their validity. 
It will turn out that these approximations have a non-negligible influence on 
the magnitude as well as the energy dependence of the pp — > ppir cross sec- 
tion. We also review the importance of contributions from pion rescattering. 
We believe that this is necessary because in the work of Hernandez and Oset 
phenomenological prescriptions for the off-shell extrapolation of the ttN am- 
plitude are applied. In the present paper a realistic meson-theoretical model 
for the ttN interaction [10] is employed. Furthermore these authors use some 
simplifications in the kinematics in evaluating the rescattering diagram which 
are avoided in the present work. We confirm qualitatively the results reported 
in ref. [8] but we find rather smaller rescattering contributions. 

Our calculations are carried out in momentum space. This allows us to retain 
all non-localities of the NN potential as well as of the ttN interaction. Fur- 
thermore we can avoid approximations at the pion production vertex. On the 
other hand we have to deal with the difficulty of incorporating the Coulomb 
interaction in the momentum-space calculations. It has been shown in the 
past [II] that the strong (on-shell) T-matrix in the presence of the Coulomb 
interaction, T cs , can be calculated reliably with a method proposed by Vin- 
cent and Phatak [12]. However, T cs is needed now half-off-shell where this 
method is not applicable. In principle one could solve the scattering equation 
in Coulomb-wave representation [13] which is however a delicate numerical 
problem. Therefore we use the following prescription for extending T cs off- 
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where k and q denote the on-shell and off-shell momentum, respectively. 
T s (k } k) and T s (q } k) are the on-shell and half-off-shell T-matrices for the 
strong interaction alone. The Coulomb penetration factor is given by 
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with M being the nucleon mass and a the fine-structure constant. We have 
checked the reliability of this prescription by applying it to the Reid soft-core 
potential [14], where the half-off-shell T-matrix including Coulomb interaction 
can be also obtained from solving the (inhomogeneous) Schrodinger equation 
[15]. The deviations occurring in the resulting pp — > ppir cross sections are 
of the order of 1%. 

We consider only the lowest partial waves in the outgoing channel; the pp pair 
is taken to be in the 1 S'o and the pion in an s-wave relative to the nucleon 
pair. From angular momentum and parity conservation it then follows that the 
initial pp state has to be 3 P - Distortions in the initial and final NN states are 
taken into account. The distorted waves are expressed in terms of the (half- 
off-shell) T-matrices in the standard way and the occurring loop integrals are 
then computed numerically. 



2 Direct pion production 

In this section we study the direct pion production mechanism shown in Fig. 
la. Here the Reid soft-core potential is employed for the distortions in the 
initial and final pp states. For this interaction model we can compare our 
results with the ones obtained by Meyer et al. [2]. The Coulomb interaction 
is neglected in this section. 

Following previous investigations we use pseudo- vector coupling for the ttNN 
vertex. This leads to the following structure for the pion production vertex 
(omitting an additional term cubic in the momenta [16] which is negligible 
near threshold), 



where p (p 1 ) is the incoming (outgoing) nucleon momentum, and e p = E p + M 



with the nucleon energy E p = \/M 2 + p 2 . u q = y m\ + q 2 is the energy of the 
pion with momentum q = p — p' . 
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In previous calculations various approximations have been made. For instance 
the a ■ (p — p') term is omitted 1 , the energies E p , £y , and u q are replaced by 
the respective masses and simplifications are made in the kinematics of the 
pion. We want to examine the effects of these approximations on the cross 
section for pp — > ppir step by step. We start with the full vertex (eq.(3)) 
which leads to the results given by the solid line in Fig. 2. 

For the non-relativistic approximation where in eq. (3) E p , E v > and u q are 
replaced by M and m^, respectively, the dotted curve is obtained. Clearly this 
leads to a 30% increase in the cross section (cf. Fig. 2). Another approximation 
commonly used concerns the evaluation of the momentum of the produced 
pion. It is done in the final NN center- of- momentum frame rather than in the 
ttNN system. This means that in the kinematical relations the reduced mass of 
the pion relative to the NN system (^(nn)tt = 2M+m ) 18 re P iace d by the pion 
mass m v . The result is an increase in the allowed maximum pion momentum 
Qmax °f arj out 7 %, which in turn causes an additional enhancement of the total 
cross section (cf. the difference of the dashed and the dotted curve in Fig. 2). 
Evidently, the two approximations discussed so far lead to an overestimate of 
the cross section by almost a factor 2. 

Finally, in most previous calculations the first term in eq. (3) has been omit- 
ted. This term is proportional to the pion momentum q and therefore it has 
been argued that it should be suppressed in the threshold region where the new 
IUCF data have been taken. However, our results shown in Fig. 2 demonstrate 
that this is true only extremely close to threshold — cf. the dot-dashed and the 
dashed curves. Already for incident energies corresponding to rj (= q^ax/ 171 -*) 
values as small as 0.2 the first term has a noticeable influence on the cross sec- 
tion and its influence increases dramatically with increasing energy. Evidently, 
the opposite signs of the two terms in eq. (3) cause a destructive interference 
leading to a decrease in the cross section. 

Fig. 2 contains also the result of a calculation by Meyer et al. [2] which is 
based on all the approximations discussed. It can be seen that our momentum 
space calculation agrees very well with their configuration space result once 
we have incorporated those approximations. 



1 This applies only to the direct production diagram. For the pion production 
vertex in the rescattering diagram this term is kept whereas the second term in eq. 
(3) is usually omitted. 
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3 Pion rescattering 



The second pion production mechanism we want to study is pion rescattering 
(Fig. lb). In most calculations the ttN scattering amplitude has been derived 
from a phenomenological effective Hamiltonian [5] 

H = 4tt— • + 47r^vl/f* • </> x <9 </> (4) 

m v m; 



where Ai, A 2 are fixed by the (empirical) Su and S31 pion nucleon scattering 
lengths [6]. The isovector term (proportional to A 2 ) does not contribute to the 
rescattering diagram in fig. lb. As Ai is very small (Ai = 0.005 according to 
Holder et al. [17]; Ai = —0.0013 following Arndt et al. [18]) it has usually 
been found that the rescattering contribution to the pp — > ppir cross section 
obtained from the effective Hamiltonian in eq. (4) is small. 

However, it is well known that the smallness of the isoscalar ttN on-shell 
amplitude is due to a strong cancellation between different isospin amplitudes. 
The situation is rather different once one moves off-shell — and it is indeed 
the off-shell ttN amplitude which enters into the rescattering diagram (Fig. 
lb). This fact is exploited in a recent investigation by Hernandez and Oset 
[8]. Using two different phenomenological off-shell extrapolations these authors 
demonstrate that the pion-rescattering mechanism can enhance the pp — > ppir 
cross section by a factor of 5 or more. 

In the following we re-examine the findings of Hernandez and Oset. In our in- 
vestigation however, a microscopic meson-exchange model of ttN interaction 
is implemented which has recently been constructed in Jiilich. This model 
includes the conventional (direct and crossed) pole diagrams involving the nu- 
cleon and the A-isobar; the meson exchanges in the scalar (<r) and vector (p) 
channels are derived from correlated 2tt exchange. Details can be found in ref. 
[10]. Here we use a slightly modified version where the form factors are energy- 
independent and the antibaryon contributions have been left out. These mod- 
ifications are made to allow an extrapolation of the model to negative energies 
as required in the present three-particle context. After readjustment of its free 
parameters this model yields a good description of low-energy ttN scattering, 
comparable to the results shown in ref. [10]. The resulting s-wave scatter- 
ing lengths are a\ = 0.173m" 1 and a 3 = —0.084m" 1 leading to a value of 
Ai = —0.001 which is in line with the value given in ref. [7]. The half-off-shell 
T-matrix produced by this model is shown in Fig. 3 at ttN threshold. 

Besides of using a microscopic ttN model there are also some differences in 
the evaluation of the rescattering contribution. In Ref. [8] only the first term 
of the pion production vertex in eq. (3) is used whereas we use again the full 
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vertex. Furthermore the four-momentum of the exchanged pion, q, is fixed 
by (f = 70 MeV and | q |= 370 MeV/c, corresponding to the particular 
kinematical situation at the pion production threshold when there are no 
distortions in the initial and final pp states. As a consequence only the half- 
off-mass-shell isoscalar ttN amplitude enters into the calculation of Hernandez 
and Oset (cf. Ref. [8]). We allow full freedom in the momentum transferred by 
the pion as it appears in the loop integrations when the initial- and final-state 
interactions are taken into account. In this case the ttN T-matrix is needed 
fully off-shell. It depends on the relative momenta of the ttN states and on 
the energy available in the ttN subsystem. 

The contribution of the rescattering mechanism to the pp — > ppir cross 
section is illustrated in Fig. 4. Note that now the Bonn potential OBEPF [19] 
is used for the initial and final pp interaction throughout and the Coulomb 
interaction is taken into account. For the pion production vertex we use the 
same parameters as for the ttNN vertex in the NN potential, namely the 
value / 2 /47r = 0.0795 for the pion-nucleon coupling constant and a monopole 
form factor with cutoff mass = 1.75 GeV. It can be seen that rescattering 
increases the cross section by a factor of around 3.5 compared to the rate 
for direct production. Our result therefore confirms qualitatively the findings 
reported in ref. [8]. However, the effect of rescattering with our microscopic 
model for the ttN T-matrix is obviously smaller. This can be understood 
if we contrast the off-shell behavior of the ttN T-matrices used in the two 
calculations, as is done in Fig. 3. We want to point out that in our calculation 
the rescattering contribution enhances the cross section even though Ai is 
negative; the effective Hamiltonian of eq. (4) leads to a decrease in the cross 
section for negative Ai [7]. 

As can be seen in Fig. 4 our predictions underestimate the empirical data 
[2,20] by a factor of 3.7. Scaling the results by this factor (dot-dashed curve in 
Fig. 4) reveals that also the energy dependence of the data is not reproduced. 
The latter feature is in contrast to essentially all previous calculations. It is 
caused by the additional ^-dependence introduced by the first term of the 
pion production vertex (eq. (3)). Thus the shape of the total cross section as 
a function of the energy is no longer determined by phase-space factors and 
the final-state interaction alone. For the sake of clarity we have repeated the 
calculation using only the second term of eq. (3) for the direct production 
and only the first term of the pion-production vertex for the rescattering 
contribution, as it is usually done in the older calculations. The result, scaled 
by a factor of 1.8, is also shown in Fig. 4 (dotted curve). 

Coming back to the direct production diagram oncemore we want to note that 
Niskanen also kept the first term in eq. (3) in his investigation [4]. He observes 
a similar interference effect as we do but in his case this interference occurs at 
somewhat higher rj values so that he can still describe the energy dependence 
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of the IUCF data. We believe that this is due to the following: Niskanen uses 
the phenomenological Hamiltonian of eq. (4) but he allows a dependence of 
Ai on the (on-shell) ttN momentum. Since Ai increases rather rapidly with in- 
creasing energy (cf. Fig. 3 in Ref. [21]) the s-wave ttN rescattering contribution 
compensates, to a large extent, the decrease in the cross section introduced 
by the first term of the direct pion production in Niskanen's work [21]. On 
the other hand, in our calculation based on the full off-shell dependence of the 
ttN amplitude the main contributions from rescattering come from off-shell 
momenta of typically 200-400 MeV/c. For these off-shell momenta the ttN 
isoscalar s-wave amplitude depends only weakly on the energy and therefore 
no such compensation occurs. 



4 Discussion and conclusion 

Our investigations reveal that several of the approximations applied in pre- 
vious studies of the reaction pp — > ppir close to pion production threshold 
are not justified. The correct treatment of the pion production vertex re- 
duces the 7T° production rate by a factor of 2. However, even more serious is 
the modification in the energy dependence of the pp — > ppir cross section. 
Whereas the calculations using the approximated vertex can reproduce the 
energy dependence of the IUCF data, the result based on the full vertex starts 
to deviate already at rj ~ 0.3 (cf. Fig. 4). This cannot be compensated by, for 
example, taking into account p-waves in the NN- and 7riV-sector since their 
contributions do not set in before values of rj around 0.5 (as we found in test 
calculations). Therefore it seems that a genuinely new mechanism is needed. A 
possible solution could be corrections from heavy-meson exchange. Indeed the 
figures shown in refs. [6,7] suggest that the contributions from heavy-meson 
exchange modify the rj dependence of the cross section to the right direction. 
However, it is unlikely that these modifications will be sufficient to remedy the 
discrepancy observed in our calculation. Another possibility which was been 
advocated in the literature [22,23] is the coupling between the final ppir , npir + , 
and o?7T + channels. In this context we would like to remark that the threshold 
of the npir + channel lies close to rj = 0.3. 

The rescattering diagram gives rise to an appreciable enhancement of the 
pp — > ppir cross section provided that the off-shell behavior of the irN am- 
plitude is properly accounted for. With our microscopic ttN model we obtain 
an enhancement factor of 3.5 which more than compensates the reduction 
which we get from the correct treatment of the vertex in the direct produc- 
tion diagram. However, we are still far away from the data. Moreover, this 
compensation is true only up to rj ~ 0.3. For larger rj values the correct oper- 
ator structure of the pion production vertex introduces an energy dependence 
which is not completely compensated by the rescattering contribution. 
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The off-shell behavior of the ttN amplitude used in the present calculation is 
completely determined by the dynamical ingredients of the ttN model and the 
fit to the ttN data. However, in the original models presented in ref. [10], which 
differs in its dynamical ingredients (inclusion of antibaryon contributions), the 
maximum of the ttN amplitude (in Fig. 3) would be roughly twice as high 
and accordingly a larger contribution from the rescattering mechanism can be 
expected. We have left out those antibaryon contributions for consistency with 
the NN model used and to permit an extrapolation to negative energies. Given 
the uncertainties in the off-shell properties of the ttN amplitude it is difficult 
at present to obtain a reliable estimate of the contribution from rescattering. 
Clearly a consistent description is needed where the NN and the ttN systems 
are treated on the same footing if one wants to make real progress. We will 
continue working in this direction. 
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Fig. 1. Mechanisms for pion production: (a) direct production, (b) pion rescatterin; 
(c) production incorporating heavy-meson exchange. 



10 



1.00 



_Q 



I 






1 i 






























X 






x \ 






/ 

X / 

// 


/ 

/ 






// ,-• 








//V 








// / / 
/ / / 








/ / / 








/ / / 






- 


; /// 








/,/ / 








- / / / 








/ .■' / 
' / 








7/7 








'.// 








i 






- 



0.1 0.2 0.3 n 0.4 0.5 

= q max /n \ 



Fig. 2. Effects of different approximations at the pion production vertex on the 
pp — ► pp7r° cross section. The Reid soft-core potential is used for the initial- and 
final-state interaction of the two protons. The solid line is the result for the full 
vertex (eq.(3)). The dotted curve is obtained when the energies of the particles are 
replaced by their masses. For the dashed line in addition the reduced pion mass 
is replaced by m^. Finally, omitting the first term in eq. (3) yields the dot-dashed 
curve. This approximation corresponds to the results obtained by Meyer et al. [2], 
which are indicated with crosses. 
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Fig. 3. Half-off-shell behavior of the ttN T-matrix at ttN threshold. The solid line 
denotes the model used in the present calculation. The dashed curve is the model 
by Hamilton [24] used in Ref. [8]. 
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Fig. 4. Results for the reaction pp — ► ppir using the model OBEPF [19] for the 
initial- and final-state interaction of the two protons. The solid curve represents our 
full calculation whereas the dashed line is obtained when only the direct production 
is taken into account. The dot-dashed curve is the full result scaled by a factor of 
3.7. The dotted curve is the full calculation based on the approximative treatment 
of the pion production vertex (cf. text). It is scaled by a factor of 1.8. The data 
points are taken from refs. [2] (dots) and [20] (squares). 
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